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ABSTRACT 

It is now generally agreed that some process prevents the diffuse gas in galaxy clusters from cooling 
significantly, although there is less agreement about the nature of this process. I suggest that cluster gas 
may be heated by a natural extension of the mechanism establishing the Mbh — o and Mbh — Mbuige 
relations in galaxies, namely outflows resulting from super-Eddington accretion on to the galaxy's central 
black hole. The black holes in cD galaxies are sporadically fed at unusually high Eddington ratios. These 
are triggered as the cluster gas tries to cool, but rapidly quenched by the resulting shock heating. This 
mechanism is close to the optimum efficiency for using accretion energy to reheat cluster gas, and 
probably more effective than 'radio mode' heating by jets for example. The excess energy is radiated in 
active phases of the cD galaxy nucleus, probably highly anisotropically. 

Subject headings: accretion - black hole physics - cooling flows - galaxies: clusters 
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1. INTRODUCTION 

Clusters of galaxies are the largest gravitationally bound 
objects in the universe. Assuming rough virial equilib- 
rium between the component dark matter, gas and galax- 
ies, with velocity dispersion a c ~ 1000 kms -1 , shows that 
within a core radius i? C orc ~ 150 kpc about the central 
cD galaxy the intercluster gas (total mass ~ 1O 14 M ) has 
a free-free cooling timescale shorter than the age of the 
Universe (see eq [3] below) . However it is by now well es- 
tablished that there is no significant mass of cooling gas 
within i?, corc flowing towards the cD galaxy, implying that 
some mechanism supplies energy to heat this gas. The 
most likely source of this energy is fairly clear: the cD 
galaxy is very massive (M c d > 10 12 M Q ) and thus probably 
has a central black hole of high mass Mbh ^ 1O 9 M0. The 
total luminous accretion energy eMc 2 released in building 
up this hole mass (here e ~ 0.1 is the radiative efficiency) 
considerably exceeds that needed to resupply the energy 
lost in radiation by the gas within i? corc (see Section 2 
below). 

However the means of transporting a suitable fraction 
of this energy to the radiating gas is far less clear. Sev- 
eral methods have been proposed, including sound waves 
and thermal conduction (e.g. Graham et al., 2008; Con- 
roy & Ostriker, 2008) and mechanical heating by jets (e.g. 
Briiggen & Kaiser, 2002; Omma et al., 2004). However 
it is uncertain how effective these mechanisms are in cou- 
pling to the gas, particularly in view of the fact that the 
cD galaxy is not active for most of the time. 

A possible alternative heating mechanism comes from 
the relation between central black hole mass Mbh and 
galaxy velocity dispersion a g observed in nearby galax- 
ies (Ferrarese & Merritt, 2000; Gebhardt et al., 2000; 
Tremaine et al., 2002), which has the form Mbh <x g\- 
Although alternative explanations exist, a promising line 
(e.g. Silk & Rees, 1998; King, 2003, 2005) suggests that 
this relation is a consequence of super-Eddington accre- 
tion on to the central black hole. This drives an outflow 
which communicates the hole's presence to the interstel- 
lar gas. At modest Eddington factors ttie = M out /MEdd 



the radiation field couples to the outflow via the single- 
scattering limit. This imparts momentum L-^dd/c to it 
(King & Pounds, 2003; King, 2003, 2005), and simulta- 
neously Compton-cools the reverse shock as the outflow 
sweeps up the galaxy's interstellar gas. Thus a forward 
shock is driven into the ambient gas purely by the momen- 
tum of the super-Eddington outflow ('momentum-driven) 
with no extra contribution from the kinetic energy ran- 
domised in the reverse shock (this would be an 'energy 
driven' flow). The dynamics of this two-shock structure 
now fix the relation between the black hole mass and the 
galaxy properties. If the black hole mass is below a crit- 
ical value C(jg the Eddington thrust L^dd/c is too weak 
to lift the interstellar gas against the galactic potential 
measured by a g . Thus the shock does not propagate out- 
side the hole's immediate vicinity and accretion can con- 
tinue. However once Mbh reaches the critical value Ca g , 
the shock attains the escape velocity ~ a g and expands to 
large radii, preventing further growth in Mbh- 

Remarkably, this simple idea gives not only the observed 
proportionality Mbh = Ca g , but also the quantitatively 
correct coefficient C = f g K/irG 2 , where f g — 0.16 is the 
gas fraction, k the electron-scattering cross-section, and 
G the gravitational constant (King, 2003, 2005). Further, 
at sufficiently large distances R c from the hole, Comp- 
ton cooling is ineffective and the extra injection of ther- 
malized kinetic energy accelerates the two shocks above 
the escape velocity cr g , driving away the remaining gas 
and fixing the baryonic mass Mbuige of the galaxy bulge. 
In the limit of modest Eddington factors (which implies 
wind outflow speeds v approaching c) one finds the re- 
lation Mbuige ~ MBii(m p /m e ) 2 a g /c ~ 10 3 Mbh between 
bulge and black hole mass, very close to observation. 

These results show that the outflows driven by super- 
Eddington accretion are very effective in communicating 
the hole's presence to the gas in the galaxy. Moreover if 
the outer (snowplow) shock reaches a large distance R c 
from the hole it strongly heats this gas because the shock 
velocity exceeds the local velocity dispersion. Here we see 
a possible connection with the cluster gas cooling problem. 
These features of super-Eddington outflows in galaxies 
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are obviously also desirable ingredients for any mechanism 
which might heat cluster gas. If the accretion energy of 
the central black hole could somehow drive a shell into the 
cluster gas, it could also reach a radius where Compton 
cooling of the reverse shock is ineffective. At this point 
the shock velocities would increase because some of the 
outflow kinetic energy is converted to heat and hence ex- 
erts pressure. This higher shock speed would exceed the 
local velocity dispersion in the cluster gas, heating the gas 
above the virial temperature and thus supplying heat as 
well as kinetic energy to the cluster gas. 

However at first sight there appears to be a major dif- 
ficulty in extending this mechanism in this way. For at 
the expected mass Mbh = Ca g the Eddington thrust 
of the hole is too weak to lift cluster gas in the clus- 
ter potential measured by a c > a g out to the radius 
R c where shock heating can be effective. This would 
instead require the considerably larger black hole mass 
Ca\ = {a c /a g ) 4 M BK ~6x lO lo of ooo M , where ctiooo = 
cr c /1000 kms -1 . 

But driving simply by the Eddington thrust L^dd/c is 
a feature of mildly super-Eddington accretion only. At 
higher Eddington ratios tHe, multiple scattering raises the 
thrust above this value, allowing a hole with only the stan- 
dard Mbh — cr g mass Ca g to drive shocks into the stronger 
cluster potential and thus heat the cluster gas out to the 
core radius i? cor e (see Section 3). This process evidently 
has high efficiency in communicating accretion energy re- 
leased near the black hole to the distant cluster gas. We 
thus have a potential explanation of the cluster heating 
problem as a natural extension of the Mbh — & problem, 
provided that we can argue that the black hole in the cen- 
tral cD galaxy should have an Eddington ratio significantly 
larger than unity. Given its privileged position this seems 
inherently plausible, and I discuss this in Section 4. 

2. CLUSTER GAS 

To fix ideas, I derive here the properties of the core gas 
in a typical cluster. For simplicity I approximate this as an 
isothermal sphere characterised by the velocity dispersion 
a c . Then the gas density at radius r is 



P 2nGr 2 
and the gas mass inside radius R is 
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where i?Mpc = R/l Mpc. Assuming virial equilibrium, 
the gas temperature is T ~ 10 s <r 2 QOO K. The dominant 
cooling process is free-free emission, with cooling time 
oc T 1 ' 2 / p. This is shorter than a Hubble time i# for 
P < Pcooi — 10 -26 criooo gem -3 , i.e. within a core radius 
(using cqn[T]) 

,Va ( 3 ) 



-Rcorc — I'^Oo'lOOO kpC 



From ([2]) the mass of this core gas is 



Mr, 



1.3 x lO 13 a^ o M (4) 

To prevent significant cooling of this gas requires an 
energy input of about 1 keV per baryon, i.e. an energy 
E h ~ (1 (keV/10 3 MeV)M corc c 2 - Kr 6 Af coro c 2 . The to- 
tal gravitational binding energy released in accreting mass 



Mace on to the central black hole of the cD galaxy is 
Eacc = tM acc c 2 , with e ~ 0.1. If ?7hcat denotes the effi- 
ciency with which this energy is used to heat the cluster 
gas we see that the total accreted mass required to prevent 
cooling is 
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(5) 



(e/0.1)?7hoat (e/0.1)?7hoat 
Clearly heating by the central black hole cannot work 
if the required mass M acCj h exceeds its total mass A/bh- 
If this is ~ 10 9 M Q we need ryhcat ^0.1. The mechanism 
described below has ryhoat —0.2. I compare this with other 
forms of heating in Section 5. 

3. HEATING CLUSTER GAS 

I suggest here that cluster gas may be heated by an 
extension of the process establishing the Mbh — & g and 
Mbh — Mbuige relations in individual galaxies, involving 
super-Eddington accretion on to the central black hole. 
The resulting outflow is roughly spherical (see below) and 
sweeps up the galaxy gas into a thin shell and tries to drive 
it out against gravity. Sijacki et al. (2007) have recently 
performed a cosmological simulation with a form of me- 
chanical feedback on cluster gas, and indeed found that it 
could prevent cooling. However much of the interaction 
between the outflow and the cluster gas necessarily occurs 
on scales not accessible to current numerical simulations. 
Here I adopt a simple analytic picture in the hope of get- 
ting some physical insight into the process. 

This type of approach is described in detail in King 
(2003, 2005). The second of these papers gives the equa- 
tion of motion of the swept-up gas shell and shows that 
this clears the accreting gas away from the central black 
hole once the hole mass reaches the critical value 



M CT = (6) 

One could follow the same procedure in considering 
the effects on cluster gas, but for our purposes a simpler 
method is adequate. We note that the weight of the shell 
of swept-up cluster gas is W(R) = GM(R)[M total {R)]/ R 2 : 
here M(R) is the mass of the shell, and M to t a i(-R) is the 
total mass inside cluster radius R, including dark matter. 
Neglecting the contribution of the cD galaxy mass, which 
is small for R ~ i? corc , this is simply M(R)/f g , since M(R) 
is just the gas mass originally inside R before the passage 
of the shock. (Note that in King, 2005, the second term 
of eqn (2) should read GM(R)[M + M totai {R)]/ R 2 , and 
the correct definition of M a immediately below eqn (3) is 
(f g K/TrG 2 )a 4 .) Using © we see that the weight W(R) is 
independent of R for large R, and has the constant value 

(7) 



W = 



G 



(the shell and total mass each increase as R, so their prod- 
uct exactly cancels the inverse-square weakening of grav- 
ity). We can now decide whether the shell reaches large 
R, and so heat the cluster gas, by comparing the weight 
W with the thrust produced by the super-Eddington ac- 
creting black hole in the center of the cD galaxy (this pro- 
cedure does not give the time dependence of the motion, 
which requires one to solve the shell's equation of motion 
taking account of its increasing inertia, cf King, 2005). 

In the single-scattering limit expected for modest 
Eddington ratios tHe this thrust is simply L^dd/c = 
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AwGM/k. This gives the expected result that the shell 
would reach large R, and thus heat the cluster gas, if and 
only if the black hole mass exceeded the value ^ with 
a c in place of a g , which as we have seen in Section 1 is 
unrealistically large (~ 6 x 10 10 M Q ). We would expect 
instead that as in other galaxies, the hole would have only 
reached the smaller value © given by the galaxy's internal 
velocity dispersion a g . 

Now let us consider the effects of an accretion episode 
with higher Eddington ratio. Shakura & Sunyaev (1973) 
consider disc accretion in this case. Their theory appears 
to apply well to X-ray binary systems, which can have far 
higher Eddington ratios than supermassive black holes in 
galaxy centers (see eq [12] below). For the well-known sys- 
tem SS433, which has to e ~ 5000, Begelman et al. (2006) 
and Poutanen et al. (2007) show that the features antic- 
ipated by Shakura & Sunyaev appear. In particular the 
total accretion luminosity is ~ ^Edd[l + In toe], and is al- 
most entirely channelled by scattering into a narrow pair 
of funnels around the disc axis, so that the outflow is es- 
sentially spherical apart from these two funnels. These two 
results suggest that highly super-Eddington accretion on 
to stellar-mass compact objects offers a plausible explana- 
tion for most if not all ultraluminous X-ray sources (ULXs: 
cf King et al., 2001; King, 2009). Most importantly for 
our purposes, most of the super-Eddington mass inflow is 
blown away from a radius i? c i rc — 9mE-Ri n /4 (where R m is 
the inner disc radius near the black hole) with mechanical 
luminosity 

1 



-Mv* ~ L E dd 
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This resulting relation v = (2L^&/ M) 1 / 2 allows us to 
estimate the thrust exerted by the accreting hole on its 
surroundings, i.e. 



Mv = (2ML E dd 
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This exceeds the single-scattering estimate by the factor 
~ (2toe/£) 1,/2 . This is potentially a lower limit to the in- 
crease, as there is a thermal pressure contribution if the 
external shock cannot cool. However Compton cooling 
still operates on the outflow, since a luminosity ~ ^Edd 
escapes isotropically through the outflow rather than via 
the funnels along the disc axis. This also means that the 
temperature profile remains flat or decreasing radially in- 
wards in the cluster center. 

The estimate © shows that the likely black hole mass 
results in a large enough thrust to heat the cluster gas if 
a significant Eddington ratio holds for some fraction of a 
cluster dynamical time -Rcorc/oc ~ 1.5 x 10 8 yr. Equat- 
ing the value ([9]) to the weight Af g a^/G of the cluster gas 
shows that the minimum black hole mass needed to get 
the shock out to the core radius is 

1/2 , 

IgK - 4 (10) 
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Put another way, a super-Eddington accretion episode can 
successfully heat the cluster gas provided that the Edding- 
ton ratio exceeds the critical value 

• r ma 180e .i 4 , . 

TO E (cnt) = a 1000 (11) 



for a fraction of the cluster dynamical time, where eo.i = 
e/0.1 and M 9 = Af BH /IO 9 M . 

4. ACCRETION 

We have seen that a sufficiently high Eddington ratio for 
the central black hole in the cD galaxy is required to heat 
cluster gas, albeit for a relatively short timescale. The 
accretion rate specified by (TTTj) is extreme - close to the 



dynamical rate ~ f g a g /G for the host cD galaxy. If the 
black hole in the latter obeys the -Mbh — c g relation, this 
rate implies an Eddington ratio 

30 (12) 



to e 



4ct„ 



where I have taken e = 0.1, a g = 300 kms -1 at the last 
step. This is in one sense reassuring, as it shows that most 
supermassive black holes in galaxy centers do not experi- 
ence very high Eddington ratios, justifying the use of the 
single-scattering limit in deriving the Mbh — o relation for 
them. Conversely, if cluster gas is heated by the process 
discussed here, the central regions of the cD galaxy must 
experience gas inflow rates of order lO 3 M yr _1 , which 
come close to destabilizing them, at least for a short time. 

There is an obvious candidate for this very rapid ac- 
cretion - the onset of the cooling catastrophe. If noth- 
ing acted to reheat the cluster gas, this would eventu- 
ally begin to flow in towards the cD galaxy at rates 

~ M COIe /t H ~ lOV^ M yr _1 . The cD galaxy must 
react long before such rates are reached. Its central black 
hole drives a snowplow shock out through the cluster gas, 
reaching the core radius and reheating the enclosed gas to 
the virial temperature in a dynamical time ~ 10 8 yr. This 
stabilizes the cluster gas and stops further infall. The duty 
cycle of the cooling-infall phase is thus of order 10 -2 . Only 
a small amount of cluster can cool before being reheated, 
so clusters typically appear to be in virial equilibrium. 

5. ENERGY BUDGET 

During an active phase of the type described above, the 
central black hole gains mass at about its Eddington rate 
for some 10 s yr, and thus typically grows by ~ 10 9 Mq. In 
return it reheats ~ 1O 13 M of cluster gas. Comparing with 
the estimate ([5]), this process uses about 5-10 times more 
than the minimum possible accretion fraction of 10~ 5 . The 
reason for this is that the central black hole manages to 
radiate about 4 times the energy it puts into mechani- 
cal luminosity (respectively LEdd[l + In toe] ~ 4iEdd, and 
^Edd) and more energy is lost in cooling. The cD galaxy 
gains ~ lO n M during an active phase, small compared 
with its current mass. We can thus write 

E q ~ 0.1ec 2 M g (13) 
for the heat input into the cluster in this form of the 
'quasar' mode from accreting mass M q on to the black 
hole. Discussions of cluster heating (e.g. Sijacki et al., 
2007) contrast the quasar mode with the 'radio' mode. 
This is motivated by observations (e.g. Birzan et al., 2004, 
Rafferty et al., 2006) which suggest that radio-loud FR I 
sources can inflate X-ray cavities. The accretion of gas 
mass M r is assumed to produce jet kinetic energy 

Ej = i] r c 2 M r . (14) 
If these jets convert their energy into cluster heating with 
efficiency rjj we get heat input 

■q r rjjC 2 M r (15) 
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into the cluster gas. Hence 

Er ^ rj r r]jM r 

E q ~ O.UM q ' 1 ' 

However several observational surveys put strict limits on 
the ratio of total jet to radiative output by AGN, or equiv- 
alently rj r M r /eM q . For example Cattaneo & Best (2009; 
see also Merloni & Heinz, 2008) find this ratio is < 0.1. 
Accordingly we find from (|16p that 

f 1 < % < 1. (17) 

Eg 

This suggests that outflows of the type considered here 
likely to be more effective than jets in heating cluster 
gas, i.e. require less black hole mass growth to produce 
the same heating effect. Thus for heating by an outflow 
the total increase AMbh in the mass of the central black 
hole is controlled by the rate at which cluster gas cools, 
i.e. AMbh - 10 -4 M coo i, so if M coo i ~ A/ C orc we ex- 
pect AA/bh ~ 1O 9 M . In principle M coo i might exceed 
M core if the reheated gas cools more quickly than before, 
i.e. in less than a Hubble time. This would then require 
multiple heating events, and thus black hole mass growth 
Mbh > 1O 9 M0. In the simple spherically-symmetric pic- 
ture adopted here this does not happen, but this conclu- 
sion should be checked by numerical simulations allowing 
for deviations from this symmetry and thus local cooling 
instabilities. In view of eqn (|17p . if this picture requires 
excessive black hole mass growth, this problem is likely to 
be worse for radio mode heating. 

6. DISCUSSION 

I have suggested that cluster gas is heated by a natu- 
ral extension of the process establishing the Mbh — CT and 
Mbh — Afbuigo relations in galaxies. The privileged posi- 
tion of the central cD galaxy means that it is intermittently 
subject to extremely high gas inflow rates. These trigger 
highly super-Eddington accretion on to the central black 



hole, which reacts by driving a shock into the infalling gas, 
efficiently reheating it and stabilizing the cluster gas. The 
duty cycle of an active phase of this type is about 10 -2 , 
so most clusters appear to be stably in virial equilibrium. 

During an active phase the accreting central black hole 
of the cD galaxy emits ~ 4LEdd ~ 4 x 10 47 ergs -1 into a 
narrow pair of cones, and ~ ^Edd ~ 10 47 ergs -1 isotropi- 
cally. If the scaling of the beaming factor with Eddington 
ratio derived by King (2008) for ULXs holds here too, an 
observer situated in these cones would infer a still higher 
apparent luminosity ~ 0.1m E x 4LEdd ~4x 10 48 ergs -1 . 
However as the cone solid angle is only ~ 10 -2 , and the 
duty cycle of active phases is also 10 -2 , it is unsurprising 
that such luminosities are not observed. An active phase 
in a cD galaxy would be observable through the isotropic 
Eddington emission. The dense outflowing wind implies a 
large photosphere, shifting the emission into the infrared, 
and a systematic search here might prove interesting. 

Comparison of the simple treatment given here with 
observation requires care. In particular a more realistic 
cluster potential must affect the ability of the swept-up 
gas shell to escape, and thus the duration of the active 
phases and the temperature structure of the cluster gas. 
Local density perturbations and the resulting cooling will 
have similar effects. These could produce shocked bubbles 
whose cooling times are significantly shorter than a Hubble 
time, as appears to be true of at least some observed cases. 
One would then require multiple heating events, totalling 
a much larger fraction of the cluster lifetime, in order to 
stave off catastropic cooling, rather than the ~ 1% total 
heating time envisaged here. Answering these questions 
requires numerical simulation. 
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